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We present the results of searches for the standard model Higgs boson at CDF in final states with 
bottom quarks. Results are derived from the complete Tevatron Run II dataset, with a measured 
integrated luminosity of 9.5 fb^' of proton-antiproton data. The searches are performed for as- 
sumed Higgs masses between 90 and 150 GeV, for Higgs bosons produced in association with 
W or Z bosons. Employing several improved techniques, these are currently the most sensitive 
searches in the world for these processes, surpassing previous CDF results by 30% beyond what 
would be expected from the addition of new data alone. Combining the search sensitivity of 
these production modes, 95% upper confidence limits on the standard model cross section times 
branching fraction are derived, yielding an observed (expected) upper limit of 4.3 (1.8) times the 
standard model prediction for a 125 GeV Higgs boson. The significance of the data relative to the 
background-only hypothesis is 2.5 sigma. 
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1. Introduction 

Understanding the mechanism for electroweak symmetry breaking, especially by testing for 
the presence or absence of the standard model (SM) Higgs boson [[I]], has been a major goal of 
particle physics and a central part of the Fermilab Tevatron physics program. Both CDF and DO 
collaborations have performed direct searches for the standard model Higgs boson |^. The new 
searches reported here include more data, additional channels, and improved analyses techniques 
compared to previous analyses. Results are derived from the complete Tevatron Run II dataset, 
with a measured integrated luminosity of 9.5 fb^' of proton-antiproton data. The searches are 
performed for assumed Higgs masses between 90 and 150 GeV/c^. 

The SM does not predict the mass of Higgs boson, mn, but the global fit of the electroweak pre- 
cision data, including recent top-quark and W boson mass measurements from the Tevatron ^, 
constrains mh to be less than 152 GeV/c^ at the 95% conference level (CL) |^. The direct searches 
from LEP Tevatron [g], and LHC results [g, ^ set the Higgs mass between 116.6 and 119.4 
GeV/c^ or between 122.1 and 127 GeV/c^ at the 95% CL. Recendy both LHC experiments gjl^] 
observed local excesses above the background expectations for a Higgs boson mass of approxi- 
mately 125 GeV/c^. Much of the power of the LHC searches comes from gg ^ H production and 
Higgs boson decays to 77, VK+W^, and Z^Z , which probe the couplings of the Higgs boson to 
other bosons. In the allowed mass range, the Tevatron experiments are particularly sensitive to 
the associated production of the Higgs boson with a weak vector boson in the bb channel, which 
probes the coupling of the Higgs boson to b quarks. 

The Tevatron collider produces proton and anti -proton collision at the center mass of 1.96 TeV 
with a record luminosity of 4.3 10^^ cm^^s^^ The Tevatron deUvered close to 12 fb^' to each 
experiment before the shutdown on 30 September 201 1. The CDF detector is a general-purpose 
detector, which provides excellent tracking, lepton identification, jets finding, missing transverse 



energy {^t) detection, and efficient multilevel triggers. The details can be found elsewhere [11]. 



2. Search Strategies 

The dominant Higgs production processes at the Tevatron are the gluon-gluon fusion (gg —?•//) 



and the associated production with aW or Z boson [|12|]. The cross section for WH production is 
twice that of ZH and is about a factor of 10 smaller than gg — )> H. The Higgs boson decay branching 
fraction is dominated by H ^ bb for the low-mass Higgs (mh < 135 GeV/c2) and by H ^ W+W' 
or ZZ* for the high-mass Higgs (mh > 135 GeV/c^). A search for a low-mass Higgs boson in 
the gg ^ H ^ bb channel is extremely challenging because the bb QCD production rate is many 
orders of magnitude lai^ger than the Higgs boson production rate. Requiring the leptonic decay of 
the associated W or Z boson greatly improves the expected signal over background ratio in these 
channels. As a result, the Higgs associated production with H ^ bb is the most promising channel 
for the low-mass Higgs boson searches at the Tevatron. 

The search strategies employed at CDF have been evolving constantly over time. We first 
maximize the signal acceptance by using efficient triggers, excellent lepton identifications, and 
powerful Z7-tagging, which can improve the signal to background ratio up to the 1% level. Then 
we use multivariate analysis (MVA) to exploit the kinematic differences between signal and back- 
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ground, which can further enhance the signal to background ratio up to the 10% level in the high 
score regions. The same strategies have been used to help discover the single-top and diboson 
processes at the Tevatron, and provide a solid ground for isolating a small signal out of a lai^ge 
background. 

For the H ^ bb signatures we look for a bb mass resonance produced in association with a 
W or Z boson where W decays into Iv or Z decays into or vv. The WH — > Ivbb is the most 
sensitive channel that gives one high Pt lepton, large ^t, and two Zj-jets. Before Z^-tagging, the 
sample is predominated by the W + light-flavor jets, which provides an ideal control data to test the 
background modeling. 

Since there are two Zj-quark jets from the low-mass Higgs decay, improving Z^-tagging is cru- 
cial. CDF use multivariate Z^-taggers to exploit the decay of long-lived B hadron as displaced 
tracks/vertices. The typical Zj-tag efficiency is about 40-70% with a mistag rate of 1-5% per jet. 



Recently CDF combined their existing Z^-taggers into a Higgs Optimized Z>-tagger (HOBIT) [ |13| ] 
using a neural network tagging algorithm, based on sets of kinematic variables sensitive to dis- 
placed decay vertices and tracks within jets with large transverse impact parameters relative to the 
hard-scatter vertices. Using an operating point which gives an equivalent rate of false tags, the new 
algorithm improves upon previous Z;- tagging efficiencies by « 20%. 

3. H ^bb Searches 

We describe the seai^ches for the low-mass Higgs boson with H ^ bb at CDF. 

3.1 WH-^lvbb 

One of the golden channels for the low-mass Higgs boson search at the Tevatron is the Higgs 
produced in association with a W boson with WH Ivbb [|l^]. We select events with one isolated 
high Pt lepton (electron, muon, or isolated track), a large missing transverse energy, and two or 
three jets, of which at least one is required to be Z>-tagged as containing a weakly-decaying B 
hadron. Events with more than one isolated lepton are rejected. For the multivariate discriminant, 
we train a Bayesian neural network discriminant (BNN) in the W + two or thi^ee jets for each Higgs 
mass, sepai^ately for each lepton type, jet multiplicity, Z^-tagging category. 

We perform a direct search for an excess of events in the signal region of the final discriminant 
from each event category. Since there is no significant excess of signal events observed in the data, 
we set an upper limit at 95% CL on the Higgs production cross section times branching ratio with 
respect to the SM predictions as a function of Higgs mass. CDF set an observed (expected) upper 
Umit at 4.9(2.8) for mh = 125 GeV/c^. 

3.2 ZH l+rbb 



Another interesting channel to search for the low-mass Higgs boson is ZH — > rl^bb []15p. It 
has a clean signature, but a low event yield due to a small branching fraction of Z — > e+e^ and 
pL^jX^ . We select events with two high Pj leptons from Z decay and two or three jets. Events 
are further divided based on lepton type, jet multiplicity, and the number of Z^-tagged jets, similar 
to WH — )• Ivbb. To increase signal acceptance we use neural networks to select loose dielectron 
and dimuon candidates. The jet energies are corrected for the missing Ej using a neural network 



3 



Direct searches for SM Higgs boson produced in association with a vector boson 



Weiming Yao 



approach. We utilize a multi-layer discriminant based on neural networks where separate discrim- 
inant functions are used to define four separate regions of the final discriminant function. There 
seems to be an excess of events in the high score signal region, but nothing statistically signifi- 
cant yet. We set an observed (expected) upper limit at 95% CL on the Higgs cross section times 
branching ratio over the standard model prediction at 7.2 (3.6) for the Higgs mass at 125 GeV/c^. 

3.3 WH,ZH ^$Tbb 

We also seai^ch for the Higgs boson in the ZH and WH channels where the Z boson decays 
into two neutrinos or the lepton from the W decay is undetected [p^. It has a large signal rate as 
well as a lai^ge QCD-multijet background. However, the final state is relatively clean, containing 
two high Ej jets and a large missing transverse energy. We require fj> 50 GeV and two Z^-tagged 
jets. We use a track-based ^^calculation as a discriminant against false If^T- In addition we utilize 
a neural network to further discriminate against the multi-jet background. The final discriminant is 
obtained for a Higgs signal by combining dijet mass, track l^j, and other kinematic variables. The 
data are consistent with the background expectations; we set an observed (expected) upper limit at 
95% CL on the Higgs cross section times branching ratio over the standard model prediction at 6.8 
(3.6) for the Higgs mass at 125 GeV/c^. 

4. Combination oiH ^bb 

We combine searches for H ^ bb in three most promising channels, and set limits with re- 
spect to nominal SM predictions [p^]. We combine the results using a combined likelihood fomied 
from a product of likelihoods for the individual channels. Systematic uncertainties are treated as 
nuisance parameters with truncated Gaussian. The common systematics between channels are the 
luminosity, detector related efficiencies, theoretical cross sections, and PDF uncertainties, which 
are treated as correlated. The instrumental backgrounds are treated as independent. Most of nui- 
sance parameters are well constrained by the data in the dominant background region and are not 
very sensitive to the initial-input values. 

To validate our background modeling and search methods, we additionally perform a search 
for Z ^ bb in association with a.W or Z boson using the same final states of the SM H ^ bb 
searches. The data sample, reconstruction, background models, uncertainties, and sub-channel 
divisions are identical to those of the SM Higgs boson search, but the discriminant functions are 
trained using the signal of Z ^ bb instead. The measured cross section of WZ + ZZ is 4.1 ± 1.3 
pb, which is consistent with the SM prediction of 4.4 ± 0.3 pb. 

To further check the consistency between data and background, we rebin the final discriminant 
from each channel in teiTns of the signal to background ratio (s/b) so that the data with similar 
s/b may be added without loss of sensitivity. The resulting data distributions after background 
subtraction are shown in figure |T] (left) and compai^ed to the expected signal of mh = 125. There 
seems to be a small excess of Higgs boson candidate events in the highest s/b bins relative to the 
background-only expectations. 

Figure ^ (right) shows the combined limit of H ^bb after combining thi^ee channels together. 
We start to exclude the Higgs mass neai^ the low-mass end between 90 < m// < 96 GeV/c^ and 
obtain an observed (expected) limit of 4.15 (1.8) at 95% CL on the SM Higgs production cross 
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Figure 1: (Left) Background-subtracted data distributions for the discriminant histograms, summed for bins 
with similar s/b, for mn = 125 GeV/c^. (Right) Observed and expected 95% CL upper limits on SM Higgs 
boson production cross section with respect to SM predictions as a function of mn- The shaded bands 
indicate the one or two sigma regions where the hmit is expected to fluctuate, in the absence of signal. The 
dot-dashes are the expected limits with the SM Higgs signal of mn = 125 GeV/c^. 

section with respect to the SM prediction for mn = 125 GeV/c^. There seems to be a small excess 
of events between 115 and 140 GeV/c^. We quantify the excess by calculating the local p-value 
for the background-only hypothesis. Figure || (left) shows the local p-value as a function of Higgs 
mass where the solid curve is for data and the dashed curve is what is expected from the SM Higgs 
production of jjih = 125 GeV/c^. At niH = 135 GeV/c^, we have a minimum local p-value=0.32% 
or 2.7 sigma, which coiTcsponds to a global p-value=2.5 sigma including a LEE factor of 2. Given 
the excess we also fit the signal cross section times the decay branching ration as a function of Higgs 
mass as shown in Figure ^ (right). We measure an associated production cross section times the 
decay branching ratio of {Owh + <^zh) 'xB{H bb) = 291+}}3 (stat-i-sys) fb, which is consistent 
with the SM prediction of 120 ± 10 fb for mn = 125 GeV/c^. 

5. Conclusion 

In conclusion, with a full dataset and many years of hard work, the CDF Collaboration has 
finally exceeded their most optimistic sensitivity projection based on 2007 summer results. The 
searches are conducted for a Higgs boson that is produced in association with a W or Z boson, has 
a mass in the range 90-150 GeV/c^, and decays into bb pair. We observe an excess of events in 
the data compared with the background predictions, which is most significant in the mass range 
between 115 and 135 GeV/c^, consistent with the Higgs-like paiticle recently observed by ATLAS 
and CMS. The most significant local excess is 2.7 standard deviations, coiTcsponding to a global 
significance of 2.5 standard deviations. This result may provide the first evidence for the Higgs 
coupling to b quarks and we are looking forward to the discovery of H ^ bb at LHC. 
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Figure 2: (Left) Background-only p-value for the combined search and the expected value assuming a SM 
Higgs signal is present as a function of m//. (Right) The best-fit cross section times branching ratio as a 
function of m//. The shaded regions show the one or two sigma bands and the SM prediction is shown as a 
smooth falling curve. The dot-dashes are what is expected for the SM Higgs signal of mn = 125 GeV/c^. 
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